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Abstract 

The gauge symmetry of the Standard Model is S U( 3) c x S (7(2)/ x U(\)y for unknown reasons. One aspect that 
can be addressed is the low dimensionality of all its subgroups. Why not much larger groups like S t/(7), or for that 
matter, S P(38) or E7? 

We observe that fermions charged under large groups acquire much bigger dynamical masses, all things being equal 
at a high e.g. GUT scale, than ordinary quarks. Should such multicharged fermions exist, they are too heavy to be 
observed today and have either decayed early on (if they couple to the rest of the Standard Model) or become reliquial 
dark matter (if they don’t). 

The result follows from strong antiscreening of the running coupling for those larger groups (with an appropri¬ 
ately small number of flavors) together with scaling properties of the Dyson-Schwinger equation for the fermion mass. 
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1. Introduction 

The Lagrangian density of the Standard Model of par¬ 
ticle physics features the gauge symmetry 

SU (3) c xSU(2 ) l xU(l)y ■ (1) 

At the hadronic scale. Quantum Chromodynamics 
(QCD, S U (3) c ) has evolved to a strongly coupled the¬ 
ory with spontaneous mass generation (and correspond¬ 
ingly, Chiral Symmetry Breaking) whereas the two 
smaller groups entail theories that remain perturbatively 
tractable, with small coupling. 

At high energies, the non-Abelian theories become 
asymptotically free and all three couplings approx¬ 
imately merge at a large Grand Unification Theory 
(GUT) scale towards which also other phenomena in 
particle physics point. 

Why these groups are symmetries of particle physics 
at collider energies is not obvious. One feature that calls 
our attention at first is the 1-2-3 succession of small 
numbers. Classical Lie groups can have arbitrary di¬ 
mensionality. Why the first three integers? It is fashion¬ 
able to resort to anthropic reasoning, perhaps within a 
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landscape of theories (“this symmetry group is compat¬ 
ible with life”), but there could also be more satisfactory 
explanations. 

In this article we adopt the view that arbitrarily 
larger symmetries could be manifest at very high en¬ 
ergy scales, but that fermions charged thereunder would 
become so massive as to be out of the reach of particle 
colliders. 

We show that if the coupling constants a s and the 
O(MeV) fermion masses are about equal for all the 
groups at the GUT scale 10 15 GeV, and compati¬ 
ble with light quarks charged under SU(3) C acquir¬ 
ing a constituent mass of about 300 MeV (so they 
are phenomenologically viable in hadron physics), then 
fermions charged under larger groups are above the 10 
TeV scale and not yet detectable. 

That is to say, fermions charged under groups of 
larger dimension than the Standard Model might exist, 
but if the coupling of those groups was similar to those 
of the SM at some GUT scale, those fermions are not 
detectable with present instrumentation. 

December 14, 2016 




We will show that the dynamical mass of those 
fermions grows exponentially with the group’s funda¬ 
mental dimension (for relatively small N c ), i.e. 

M(0V oc e Nc x <9(iV“ ticaI - Nf) (2) 

and then increases more slowly for larger N c , saturat¬ 
ing towards the GUT scale (where all are equally light 
by construction). The Heavyside step function in fla¬ 
vor limits the validity of the result to fermions whose 
flavor degeneracy is smaller than a certain critical value 
at which the vacuum polarization becomes insufficiently 
antiscreening (and beyond which dynamical chiral sym¬ 
metry breaking ceases). This is further discussed below 
in subsection 3.2. 

To establish the result shown in Eq. (2), we will find 
rescaled solutions of the mass Dyson-Schwinger equa¬ 
tion that allow us to avoid difficult numerical integration 
over large intervals of momentum. We will use these 
solutions in conjunction with a perturbative analysis of 
the highest energy scales, where a s is small. The key of 
the analysis is to note that the scale at which the cou¬ 
pling constant times the relevant color factor becomes 
sizeable, so that the DSE needs to be employed (which 
for concreteness we will take as (Cfa s ) = 0.4) is larger 
for larger groups due to the increased antiscreening in 
Yang-Mills theories, so that the fermion mass runs for 
larger intervals and thus becomes much larger at p = 0. 

In section 2 we introduce and simplify, following 
standard theory, the DSE for the fermion propagator. 
There, in subsection 2.2, we will already change the 
group under which the fermions are charged and ob¬ 
serve, numerically and at fixed cutoff, that the solu¬ 
tions for larger groups seem to be simple rescalings of 
the known 5(7(3) solution. In subsection 2.3 we will 
change to the MOM scheme to avoid the inconvenients 
of cutoff solutions. Section 3 takes us to the highest 
energies where the use of perturbation theory is appro¬ 
priate, and we will briefly recall antiscreening and per¬ 
turbative mass running in non-Abelian Yang-Mills the¬ 
ories. 

The crux of the article is then section 4, where the 
scaling properties of the rainbow-ladder DSE are com¬ 
bined with the perturbative analysis to yield our main 
result, shown in figure 11: that the fermion mass be¬ 
comes very large for larger groups, and that it scales for 
moderate N c as in Eq. (2). Further discussion spans sec¬ 
tion 5. The appendix is reserved for mathematical detail 
(computation of the group color factor Cp are reported 
there). 


2. Some properties of spontaneous mass generation 

The mass function plays a central role in gauge theo¬ 
ries coupled to fermions and their uses for phenomenol¬ 
ogy. A brief summary discussing several subtleties and 
identities is given in [1], 

We want to adopt the simplest possible Lorentz- 
invariant model that exposes the physics. The Nambu- 
Jona-Lasinio model is a practical option to demonstrate 
spontaneous mass generation, but its contact-interaction 
structure cannot be used at high energies, where the cou¬ 
pling is not transparently related to the running coupling 
of the underlying non-Abelian theory. 

Next in difficulty is the rainbow approximation to the 
Dyson-Schwinger equation [2] of the fermion propaga¬ 
tor in the gauge theory, so we settle to it [3]. While a 
very basic approximation, the simplicity of the scenario 
we propose does not require more sophisticated many- 
body methods. Rainbow-ladder approximation is still 
widely used for exploratory studies of beyond the stan¬ 
dard model physics [4]. 


2.1. Dyson-Schwinger equation for a fermion propaga¬ 
tor 

The free propagator of a fermion with current mass 
m c is denoted as 

So(p 2 ) = ■ (3) 

f - m c 

The full propagator is usually parametrized as 


S( P 2 ) = 


i 

A(p 2 )i/> - B(p 2 )' 


(4) 


but to expose spontaneous mass generation it is suffi¬ 
cient to consider a simplified ansatz with A(p 2 ) = 1 and 
running mass B(p 2 ) = M(p 2 ) = M p . 

The Dyson-Schwinger equation (DSE) for this full 
propagator. 


S-V) = S 0 V)-2(/r) (5) 


may be written down as an identity in the field the¬ 
ory, but can pedagogically be deduced as a resumma¬ 
tion of perturbation theory. The rainbow resummation 
avoids all diagrams with vertex corrections, counting 
only those of the type depicted in figure 1 . 

After standard manipulations 1 , the DSE takes the 
well-known form 


M p = m c + 


Cf r 

7T 3 JO 


q 3 dq 


\q\ 


Mq g 2 P° 
2 . ^ # 


+ M 2 


p-r 


( 6 ) 


1 Tracing over Dirac matrices, performing a Wick rotation to Eu¬ 
clidean space q° —» iq°, p° —> ip°, f d 4 q —» if d 4 qE, and employing 
4D spherical coordinates. 
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Figure 1: Resummation of the rainbow diagrams (with perturbative 
gauge boson propagator and fermion-boson vertex) leading to the 
DSE for the fermion propagator in rainbow approximation. 

where Cf is the color factor (or Casimir of the group’s 
fundamental representation) which is the object that we 
will vary in this investigation. Also seen are g, the 
fermion (non-Abelian) charge; and the Feynman-gauge 
gauge-boson, or for short even beyond QCD, “gluon” 
propagator 

- lD((p - qf)q pv = — (7) 
(p - q) 

averaged over 4-dimensional polar angle, 

1 

Vl - x 2 D p - q dx = D° p _ q . (8) 

l 

This f/VV-indcpcndcnt) gauge boson propagator is taken 
to be perturbative, though if need be, this can be cor¬ 
rected in future work to achieve better precision (see [5] 
for a very brief outline of the current estimates in non- 
Abelian gauge theory, and [6] for more extended dis¬ 
cussion). The use of the same propagator for all N c is 
supported by independent studies [7]. 

To solve the DSE we discretize the variables p, q and 
the function M, so the g-radial and x-angular integrals 
become discrete sums (needing regularization as they 
are divergent at large q), and linearize M = Mq + m 
where Mo(p 2 ) is a guess and m(p 2 ) the unknown cor¬ 
rection returning the correct solution M{p 2 ). Expand¬ 
ing Eq. (6) to first order in m provides a linear system 
for m{p 2 ) solved with a linear algebra package. The im¬ 
proved M(p 2 ) is used as a new guess Mo(p 2 ) and the 
procedure iterated until m — 0. 

2.2. Mass generation at the hadron scale (with cutoff 
regularization) 

To show the reaction of the DSE Eq. (6) to chang¬ 
ing the group, we first study the hadronic scale cutting 


off the q integral at A = 10 GeV. We take the (cutoff- 
dependent) current mass m c = m( A 2 ) = 0 for the free 
fermion to vanish, and solve for M(p 2 ) at smaller scales, 
so the entire mass function is here dynamically gener¬ 
ated breaking the global chiral symmetry. 

To be specific, in the calculations shown in figures 2 
and 3, the coupling g is taken to be the same for all 
groups and fixed by demanding that the quark mass 
for S U (3) be 300 Me V, as corresponds to the observed 
QCD quarks. This results in a value g - 15.1 with the 
cutoff fixed at 10 GeV. 

This value of g amounts to a s — 18, much larger 
than one naively expects in QCD. This is due to several 
reasons, among them having set m c = 0, which sup¬ 
presses M( 0) a moderate amount; having fixed A = 10 
GeV, which restricts the range of running mass a bit; and 
saliently, the use of a bare qqg vertex. Since chiral sym¬ 
metry breaking has to be simultaneous in all Green’s 
functions [8] and there is feedback between them, our 
use of a bare vertex underestimates the extent of chi¬ 
ral symmetry breaking, requiring a larger a s for equal 
M( 0). 

We can think of this larger g as simply the product gV 
with V a vertex strength factor. For S U( 3) this factor is 
7.7, and for other groups it scales as V(N C ) = Vsu( 3 )y, 
which is the leading N c behavior of the vertex one- 
loop corrections (specifically, the non-Abelian correc¬ 
tion). There is a vast literature on vertex corrections 
that dates back decades, see e.g. [9] or [10], so we ab¬ 
stain from further investigation here as this would carry 
us too much off topic. But it is clear that the rainbow- 
ladder approximation is just a first approximation to the 
physics. 

The results for a couple of special groups and for all 
the classical Lie groups (S U(N C ), SO (N c ) and Sp(A c ), 
with N c being even for the later) clearly show mass 
functions that seem to be rescalings of one another upon 
changing the group dimension 2 , with mass generation 
almost directly proportional to the fundamental dimen¬ 
sion of the group. 

Let us now concentrate on the “constituent” mass 
M( 0) seen at lowest energies, while varying the color 
number N c and the group families. For this we extract 
the first point of each M(p 2 ) function and plot the out¬ 
come in figure 4. 

From the figure, it stands out that for U( 1) and 5 U(2) 
there is no chiral symmetry breaking, i.e. M( 0) = 0, for 
the same coupling intensity that generates the 300 MeV 
quark mass in S U( 3). Past dedicated S U( 2) (and also 


2 we will elaborate on this property later in section 4. 
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Figure 2: Mass function for the special unitary SU(A C ) (top) and 
orthogonal SO(N c ) (bottom) with momentum integral regularized at 
A = lOGeV. In this fixed momentum interval the constituent mass 
grows nearly linearly with the group dimension (of the fundamental 
representation). 


Figure 3: Mass function for the symplectic groups Sp (N c ) (top) and a 
couple of special groups as indicated (bottom) with momentum inte¬ 
gral regularized at A = lOGeV. Again, in this fixed momentum inter¬ 
val the constituent mass grows nearly linearly with the group funda¬ 
mental dimension as in fig. 2 . 


G2) lattice studies[7, 11] found that the general struc¬ 
ture of the Green’s functions is similar to the 5(7(3) 
case, for commensurate but larger coupling (presum¬ 
ably to make up for the reduced color factors) at a low, 
hadronic scale. Our setup, and thus our result, differs 
in that the couplings are equal at a high-energy scale so 
the coupling for smaller groups is much smaller at the 
lower scale. 

A related, dedicated study [12] shows how lowering 
the antiscreening of QCD eliminates dynamical mass 
generation. 

Beyond (7(1) and S (7(2), we find no mass generation 
for G 2 (N c = 7) and f 4 (N r = 26), both with a relatively 
small color factor Cf- 1 in spite of their large dimen¬ 
sion; and also for S 0(N C ) with N c = 1 to 5 and for Sp( 2). 
For all these groups, an explicit fermion mass m just 
yields an M(p 2 ) that slightly separates from the pertur¬ 
bative value without really yielding symmetry breaking. 

For the rest of the classical groups, where symme¬ 
try breaking is apparent, the dependence of M( 0) on the 
defining dimension N c is seen to be rather linear. This, 
as we will see, happens because we have integrated over 
the same momentum interval (0, A). 


From the linear dimension of the leading divergence 
in the DSE one can also deduce that M( 0) oc A, which 
can anyway be checked numerically as shown in fig¬ 
ure 5. 

After this warmup, we have shown that mass genera¬ 
tion at the hadron scale is insufficient to expel fermions 
charged under large groups from the spectrum. This is 
no longer true when considering high-energy physics, 
where running over large momentum swaths is in¬ 
volved. But before proceeding, we note that cutoff reg¬ 
ularization is inadequate (now that the highest scale will 
be pushed to 10 15 GeV), so we first introduce an appro¬ 
priate renormalization scheme in the next subsection. 

2.3. One technical improvement: momentum subtrac¬ 
tion scheme 

There are many reasons to improve on simple cutoff 
regularization, among them preserving Lorentz invari¬ 
ance and exposing renormalizability. To characterize 
the quantized theory we need a renormalization scale 
p at which the couplings a s = g 1 /An are to be cho¬ 
sen. To achieve this, we will adapt a variation of the 
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Figure 4: Dependence of the constituent mass M( 0) with the color 
number N c under a cutoff regularization with A = lOGeV. For a given 
classical group family, the dependence is rather linear. 


Momentum Subtraction Scheme or MOM often used 
in this subfield of Dyson-Schwinger equations. Since 
we will later, in our perturbative analysis, employ only 
1 -loop running of masses and coupling constants, we 
can take the renormalization group coefficients /3 and j 
to be the same as in the more usual Modified Minimal 
Subtraction Scheme (MS), as they are equal to one loop 
(see [13]). 

The first step is to introduce adequate renormalization 
Z(A 2 ,/i 2 )-constants that absorb any infinities or, once 
regulated, any dependence on the cutoff A, 

S- 1 (p 2 ,fi 2 ) = Z 2 So\p 2 )-^(p 2 ,^ 2 ), 

£ s ig 2 C F J ( q 2 ,p 2 )y v D((p - ^) 2 ,^ 2 ) , 

namely Z 2 for the wavefunction renormalization and Z m 
for the bare quark mass. We do not calculate vertex cor¬ 
rections nor loops involving ghosts in this article since 
they are an unnecessary complication for the physics 
exposed, so we need no additional Z constants beyond 
those of the bare quark (inverse) propagator, Sq 1 (p 2 ). 
Therein, the relation between the (cutoff dependent) un¬ 
renormalized mass m c ( A 2 ) and the renormalized mass 
at the renormalization scale m R (p 2 ) is [14] 

m c (A 2 ) = Z m (AV)m«(/4 2 ). (10) 



Figure 5: Dependence of the SU( 3) constituent mass M(0) with the 
cutoff. This very nicely linear relation would get modified in a more 
sophisticated truncation of the gauge theory where the gluodynamics 
generates an additional scale (a gluon mass-like parameter that cuts 
the propagator in the infrared, another topic on which there is a large 
literature). We stay with a strictly massless gauge boson propagator 
as in Eq. (7) through the entire article. 


Should we lift the restriction A — 1, the renormaliza¬ 
tion of the wavefunction would entail Aq'(/? 2 , A 2 ) = 
Z 2 A~ l (p 2 ,p 2 )\ though while we maintain it, then also 
Z 2 — 1 and the only needed renormalization condition 
is to fix the mass at p 2 = p 2 . The DSE for the mass 
function is then formally 

M(p 2 ) = Z m m R (p 2 ) + Z M (p 2 ,p 2 ); (11) 

evaluating it at p 2 = p 2 and subtracting both, we obtain 

M(p 2 ) = M(p 2 ) + ZJp 2 ,p 2 ) - I. M (p 2 ,p 2 ), (12) 


in terms of finite quantities alone. Thus, the resulting 
MOM equation is 

8 2 C f 


M(p 2 ) = M(p 2 ) 
M(q 2 ) 


f 


q dq 


(D° - D° ) 

2 + M 2 (q 2 ) p - q ^ 


03) 


(with p parallel to p), that is, 

8 2 C f 


M(p 2 ) = M(p 2 ) ■ 


r*oo 

I q>dq L 


(lx Vl 


1 


1 \ M(q 2 ) 


) M 2 (q : 


) + \q\ 2 


r- (14) 
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If the q radial integral in this equation is cutoff at A » 
(p,p), it is easy to see that asymptotically, for p and p 
parallel, 

dM(p 2 ) M(A 2 )(p - p) 

~sK~ “-X?- (15) 

so that for large A and M growing slower than quadrat - 
ically at large momentum, M(p 2 ) stops depending on 
the cutoff, renormalization is achieved and M{p 2 ) alone 
determines the function for values of p smaller than p. 

We again fix (for all groups) g — 15.07 at p — 10 GeV 
so that for S (7(3) the constituent quark mass is M{ 0) = 
300 MeV once more. We impose the renormalization 
condition M(p 2 )suo ) = 5.7 MeV for all groups. The 
tail of the mass function for the group SU(3) approaches 
zero asymptotically, as shown in figure 6. 

Also shown are mass functions for S (7(4) and 5 (7(5) 
that are seen to change sign. This is not necessarily that 
the computer code has found the excited, sign-changing 
solutions of [15, 16, 17]. Instead, what it shows is that 
the self-energy at p subtracted in Eq. (14) is very large 
and overcomes the smaller self-energy computed at p as 
well as the smaller mass chosen at 5.7 MeV. This sim¬ 
ply reflects a renormalization point p that is too low for 
the higher groups, before the perturbative behavior sets 
in (but we want to compare the three functions at the 
same point), so we are not only subtracting the ultravio¬ 
let divergence but also large finite-/? contributions. This 
suggests, as we will soon effect, to move the renormal¬ 
ization point of the larger groups to a much higher scale 
where the coupling is weaker. 

Ignoring that sign for now, the solutions are seen to 
be similar in shape to the ones obtained with the cutoff 
method. Turning now to the deep infrared, we conclude 
that the outcome is equivalent to that obtained in subsec¬ 
tion 2.2, with M(0) scaling in proportion to N c if only 
the hadron scale is considered, so we have achieved a 
very simple renormalization that allows us to proceed 
to higher scales. 


3. Treatment of the high-energy running mass 
within perturbation theory 



Figure 6: Mass functions for some SU(N C ) groups in the MOM 
scheme with renormalization point p = lOGeV, where we have cho¬ 
sen M(p 2 ) = Msu( 3 )(p 2 ), the latter such that Msu{ 3)(0) = 300 MeV. 
To compare the very different growth of the three mass functions even 
at low scales, we have chosen them equal at a very low p so that the 
S U{ 4), S U( 5) ones eventually become negative at high energies. This 
is of course unphysical, and just means that M(p) should naturally be 
chosen larger because chiral symmetry is already broken. We never¬ 
theless find the plot instructive. 


3.1. Running coupling and mass 

In that energy regime, the running of the mass and 
coupling constants can be followed in perturbation the¬ 
ory, as long as a s remains small. Up to one loop, we will 
need the f3\ and y\ coefficients of the /3-function and of 
the anomalous mass dimension, respectively 


da ^ 

(5(a s ) = ~p— 
dp 

= Pi a 2 s + Pi a l + 

with a(s ) = —, and 


p dm 

j(a s ) = -— 

m dp 

= y\a s + yia 2 s + ... 


06) 


(17) 


We now extend our study to the Grand Unified The¬ 
ory scale at 10 l5 GeV. Several physics coincidences 
point out to some dynamics taking place at that scale, 
for example the see-saw Majorana mass scale in neu¬ 
trino physics, and most important for this work, the 
approximate coincidence of the coupling constants of 
the Standard Model gauge interactions at that scale 
(see [18] for an introductory review). 


We will, for simplicity of the argument, consider that 
there is only one fermion flavor charged under each of 
the color groups, so that we may set Nf — 1. Follow¬ 
ing [19, 20], we have 

Pi = l(ll!V c -2 N f ), (18) 

D 

n - \c F , (19) 
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and the solutions to eqs. (16) and (17) is obtained after 
integrating once. 


f 

J a s 

f 

Ja , 


a s (Pl) 

■a s (M2) 

l) 


da s >ui 

-= In —, 

P(a s ) /J2 

(20) 

y(a s ) m(p 2 ) 

aa s -= In-, 

' P(a s ) m{n i) 

(21) 


from which follow the well known forms 

1 


a s (fi 2 ) - a s (fii) 


l + S^fdgjln^’ 

n ^ 


( 22 ) 


and 


m s (ju 2 ) = m s (ji 0 


x a 

Mi 


1 + In a 

n ri n 


(23) 


In what concerns our study, it is worth remarking that 
groups of equal dimension in different families have dif¬ 
ferently running masses (for equal and low flavor num¬ 
ber, in our estimates Nf = 1). This is in spite of the run¬ 
ning of a s depending on the chosen group only through 
its defining dimension N c (of course, equal to the ad¬ 
joint Casimir Ca )■ The reason is that the actual color 
factor that appears exponentiating the fermion mass in 
Eq. (23) is the Casimir Cf of the fundamental represen¬ 
tation, which is different for two groups belonging to 
different families even if they have the same dimension 
(in short, equal N c does not imply equal C/- [/V r ]). 

These running masses are depicted in figures 7 and 8, 
that already hint at much heavy fermions even in pertur¬ 
bation theory. 

Returning to the running coupling, [f is positive for 
non-Abelian Yang-Mills theories (N c + 1), so a s (jj. 2 ) 
decreases logarithmically for fj 2 > >ui, and asymptotic 
freedom is manifest. Running in the opposite direc¬ 
tion towards lower energies, the intensity of interaction 
increases until a Landau pole A is hit (not to be con¬ 
fused with the earlier cutoff), where the denominator of 
Eq. (22) vanishes. 


A = yiqe . 


(24) 




Figure 7: Running mass M(p 2 ) as a function of N c from perturbation 
theory running at one loop from the GUT scale. Here we depict the 
classical unitary and orthogonal group families. 


Much earlier than that pole, these analytical formulae 
cease to be applicable and must be substituted by resum¬ 
mation, e.g. by DSEs. The Landau pole is of course a 
notorious feature of perturbation theory, that is avoided 
in other approaches. In Analytical Perturbation Theory, 
for example, a s saturates at low energies [21]; Dyson- 
Schwinger equations studying the gluon-ghost sector 
of Landau-gauge QCD concur [22]; and generally one 
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does expect a flattening of a s at low scales, yielding a 
conformal window [23]. 

Therefore we need to match the high-energy treat¬ 
ment, that can be handled in perturbation theory as just 
explained, with the earlier DSE treatment at some scale 
m(i u 2 ), which is the object of the next section. 



log 10 (p) (GeV) 


Figure 8: Running mass M(p 2 ) as a function of N c from perturbation 
theory running at one loop from the GUT scale. Here we depict the 
classical symplectic group family and some exceptional groups. 


3.2. Effect of the number of flavors 

The reader will have noticed that Eq. (18) depends on 
the number of flavors, which we have taken as Nf - 1 
for the numerical examples (in lattice language, this is 
the “quenched approximation”). However, as it is well 
known, if there is a sufficiently large fermion degener¬ 
acy, which in one-loop perturbation theory as encoded 
by that equation is Nf = —the vacuum polarization 
becomes screening instead of antiscreening (the sign of 
Pi changes). 

The number of flavors necessary for this screening in 
S U( 3) is 17, and for S U( 4) it is 22, and larger yet for 
higher groups, which seems a rather large degeneracy. 
However, for smaller Nf one may have an antiscreening, 
yet too weak, interaction that will fail to trigger dynami¬ 
cal chiral symmetry breaking and thus a nonperturbative 
fermion mass. 

Estimates of the critical flavor number beyond which 
chiral symmetry breaking ceases have been provided in 
the literature. Closest in spirit to our work are those 
from the DSEs [24] as well as the Renormalization 
Group Equations [25]. The DSE estimate in [24] is, for 
SU( 3), A amical = 8 + 1. The second work quotes nu¬ 
merical estimates that are compatible within the error, 

^critical = 1! + 2. 

Because of Eq. (18), it is plausible that A™ tlcal oc N c , 
so that the number of flavors necessary to overturn chiral 
symmetry breaking keeps growing (so that, for example, 
for N c - 4 we would have iV™ tlcal = 11+2 or 15 + 3 
respectively). 

The existence of this critical number of flavors jus¬ 
tifies the 6»(iV™ tlcal - Nf) factor in Eq. (2): above that 
number, M(0) becomes of order the current mass m c and 
depends only radiatively on N c . Nf acts as the parame¬ 
ter of a quantum phase transition and our results apply 
only to the broken symmetry phase. 

As a digression, for Nf close but above A^ ntical , 
our rainbow-ladder approximation in section 2 yields 
Miransky scaling (see e.g. [26]), by which M( 0) oc 

Aexp((const.)/( ^Nf - A™ tlca1 )). Beyond rainbow- 
ladder, this exponential becomes modified to a power- 
law; the window above iV™ tlcal during which these criti¬ 
cal behaviors are active is however very small [25] about 
a few percent of iV™ tlcal (see fig. 5 of that work), so that 
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for Nf = AY mical - 1 we can safely consider ourselves in 
the broken phase. 

In conclusion of this subsection, though most of 
the considerations in this article are for a flavor- 
nondegenerate fermion charged under the various Lie 
groups, they can actually be extended to Nf of modest 
size below A“ tical . 

4. Mass running from both high and low energies 

In this section we seek to combine the perturbative 
running at large scales with the DSEs at lower momenta, 
to obtain a picture which, even if crude, is global and 
allows a general statement to be produced. We start the 
perturbative renormalization group running at pcur = 
10 I5 GeV, where we fix 

£ *s(Hgut ) = 0.017 , mificur) = 1 MeV . (25) 



Figure 9: One-loop running coupling for the SU(N C ) (N c > 3, Nf = 
1) family of special unitary groups. For other families the running 
is identical since Eq. (22) depends, through /3\, only on the group 
fundamental dimension. All couplings are chosen to be identical at 
the GUT scale 10 15 GeV. 


This fermion mass is chosen to broadly reproduce the 
value of the S t/(3)-colored quark mass, that under 
isospin average, is taken [18] to be about 


m(2GeV) = 


m u (2GeV) + mj(2GeV) 
2 


3.5MeV .(26) 


As for the coupling constant, the one corresponding to 
SU(3) is precisely known at the Z-boson scale, p = 
M- =* lOOGeV (91.2GeV), where a s (M z ) =* 0.12. 
Running to one loop and with only one fermion flavor 
charged under each group (shown in figure 9) requires 
an a s at the GUT scale that is somewhat smaller than 
the usually quoted value a s (GUT ) - 0.025. But all to¬ 
gether we seem to differ by a moderately small factor 
which does not affect our main argument. 

We use the perturbative formulation encoded in 
Eq. (22) from p\ = pour down to cr = p 2 where cr rep¬ 
resents the point where perturbation theory breaks and 
non-perturbative methods are required. For S U{3), this 
point is characterized by a s = 0.3, where we decide that 
perturbation theory must break down quickly. The ac¬ 
tual combination appearing in the DSE is i,’ 2 CY o: s Cf- 
Therefore, a s = 0.3 for S U( 3) is equivalent to Cya s = 
^ x 0.3 = 0.4. From that point on, we freeze a s to a 
constant value and employ Dyson-Schwinger methods 
to treat the fermion mass. 

In fig. 10 we represent our complete approxima¬ 
tion for the quark mass function in S U{3). We match 
the perturbative and DSE solutions continuously (ob¬ 
taining a smoother matching is possible by employ¬ 
ing resummed perturbation theory on the high energy 
side [27]). 


If we now increase the dimension of the group G, 
the matching point cr where Cf-a J (cr) = 0.4 and a non- 
perturbative treatment starts to be required moves much 
to the right of the plot to higher scales, 

cr = p CUT X e^l^r^Gyr)) _ (27) 

The exponent being negative and proportional to N~ l , 
increasing N c moderately provokes an exponential in¬ 
crease in the scale. When N c becomes large, cr —> pcuT 
saturates and basically all further groups require non- 
perturbative treatment from early on. 

Integration to such large scales with an appropriate 
grid is time consuming; it can be avoided by noticing, 
for example after a glance at figures 2 and 3, that given 
a solution to the DSE’s, one can easily find rescaled so¬ 
lutions. In those figures the color factor induced the 
rescaling, but now the rescaling will rather be forced 
by cr, the point where we start numerical integration to¬ 
wards lower values of p. 

We will obtain solutions for groups of large dimen¬ 
sion from that of N c - 3 shown in figure 10. To show 
that this is possible analytically, perform a scale trans¬ 
formation 

p 2 -> '*y 

o- 2 -> A V, (28) 

on the DSE, where A is a contraction factor that will 
map the mass function of an arbitrary group to that of 
S U{3). We can always change the dummy integration 
variable q 2 —> A 2 q 2 , and the integration measure picks 
up a Jacobian d A q —» A 4 d 4 q. 
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Figure 10: Mass function for the SU(3) group obtained matching the 
numerical solution of the DSE to the perturbative one for a s (o-) = 0.3. 


With this rescaling, the DSE in Eq. (14) becomes 
8 2 Cf 


M{A 2 p 2 ) = M^o 2 ) ■ 


M(A 2 q 2 ) 


f 


A 4 q 2 dq 


D° 

p-q 


A 2 q 2 + M 2 {A 2 q 2 ) 


If 


A 2 


A 2 


(29) 


It is easy to find the modified M that satisfies this 
rescaled equation. Taking simply M(A 2 p 2 ) = AM(p 2 ) 
we indeed recover Eq. (14) so if M solves the former, 
M solves the newer, rescaled one; and the corresponding 
relation for the constituent masses is simplest, 


M{ 0) = 


M( 0) 

A 


(30) 


Figure 11: Dynamical mass M{ 0) as function of N c from 3 to 12 ob¬ 
tained by matching perturbation theory and DSE when {Cf&s) = 0-4, 
and obtaining the DSE solution by rescaling that of SU(3). Be¬ 
cause Cf (and more so Ca ) are basically proportional to N c , there 
is not much difference between a full non-Abelian theory and a quasi- 
Abelian truncation in rainbow-ladder approximation, and the scaling 
is qualitatively similar for the three families of classical Lie groups. 
This stops being true for Nf nearing the critical value, when the ex¬ 
ponential Miransky scaling (see subsection 3.2) of the quasi-Abelian 
truncation changes to a power-law.) 


is charged, the exponential in Eq. (2) for moderate N c . 
This is the reason why fermions charged under a large 
group are expelled from the low-energy spectrum, all 
things being equal at the GUT scale. 

Carrying out the rescaling for several values of N c 
leads to the dynamical mass M{ 0) dependence on N c 
depicted on figure 11. 

5. Discussion and outlook 


We put this scaling property of the rainbow DSE to 
use immediately. Taking A as the ratio of saturation 
points where a s = 0.4 /Cf, 


o group 

VSUO) 


(31) 


the mass function rescales in the same way: 

Mgroupif) _ 

Msu( 3)(0) 


(32) 


or simply put, eliminating the auxiliary A, 

MgroupiO) O' group 
3Tst/(3)(0) CTSUO) 


(33) 


This is a central result. When combined with the expo¬ 
nential growth of the saturation point in Eq. (27), we ob¬ 
tain our advertised dependence of the fermion mass with 
the fundamental dimension of the group under which it 


The combination of two methods (perturbation the¬ 
ory and the Dyson-Schwinger equations) has allowed 
us to show that fermions charged under a large group, 
if their coupling is equal to the smaller-dimension ones 
that appear in the Standard Model at the GUT scale 
10 15 GeV, are much more massive than the ones we 
see. In fact, should there exist fermions charged un¬ 
der S 1/(4) or a group of equal dimension, they would 
appear in the 10 TeV region, though we cannot pin¬ 
point them to better than order of magnitude estimate 
because of the crude approximations we have made, but 
they would not be far out of reach of mid-future ex¬ 
periments. Perhaps precise calculations in the near fu¬ 
ture can address this dimension-4 group to predict the 
mass at which S f/(4)-charged fermions appear. One 
can conceive a combination of methods coming together 
to obtain a good prediction: lattice QCD techniques 
that have already been demonstrated for groups larger 
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than in the SM [28, 29, 30, 31], scaling properties of 
full DSEs or the Exact Renormalization Group Equa¬ 
tions [32, 33, 34], and multiloop perturbation theory. 

It already appears from our simple work that groups 
yet larger might just endow fermions with a mass not 
detectable in the foreseeable future. 

Should these superheavy fermions be coupled to the 
Standard Model, they would have long decayed in the 
early universe due to the enormous phase space avail¬ 
able. Were they to exist and be decoupled from the SM, 
they would just appear to be some form of dark matter. 

In addressing the spectrum of Beyond-SM theories 
one can worry that spontaneous mass generation may 
break any extant global chiral symmetries and give rise 
to presumably unseen Goldstone bosons equivalent to 
QCD’s pions. To dispel doubts, let us recall the Gell- 
Mann-Oakes-Renner relation [35] 

M\fl = -2 m q (qq) (34) 

relating quasi-Goldstone boson mass and decay con¬ 
stant to fermion mass and condensate. The dependence 
with the typical scale of symmetry breaking is f n ~ A, 

(qq) -A 3 , and therefore M n ~ Am q (A ) (note that 

m q (ncuT ) = O(MeV), and it is bigger at A). This 
puts the pseudo-Goldstone bosons out of reach of con¬ 
temporary experiments, except perhaps for the group 
S 1/(4) and equal-dimension ones. In detailed modeling 
one can also try to arrange for quantum anomalies lift¬ 
ing the necessity of unwanted Goldstone bosons, such 
as QCD’s q'. We abstain from attempting this at the 
present time. 

We have shown that the fermion mass for groups 
slightly larger than S 1/(3) grows exponentially with N c , 
because the mass satisfies the same scaling relation than 
the saturation point, <x, of the coupling constant a s 
(which is obviously a proxy for some more sophisti¬ 
cated saturation mechanism), and this point grows ex¬ 
ponentially with N c according to Eq. (27). 

In our discussion there is a degree of arbitrariness: we 
have assumed that the coupling corresponding to larger 
groups at the GUT scale, which is totally unknown, is 
the same for all groups (after all, that is the meaning 
of GUT). If this hypothesis is lifted, one can of course 
find arbitrary results. Just like QED with stronger 
coupling can generate mass spontaneously [36], very 
large groups with sufficiently small coupling at the 
GUT scale 3 would not generate it and we would have 
fermions charged under the oddest groups at current col¬ 
lider scales (which does not seem to be the case). We 


also emphasize that our discussion has focused on one 
or at most few new flavors. If a large N c is accompanied 
by a very large Nf one can overcome the gauge-boson 
antiscreening with fermion screening. Our conclusions 
then need to be revised. 

We insist once more that the couplings for all groups 
are taken to be the same at pcur, and we do not suppress 
them as in t’Hooft’s counting [37] with g oc which 
may induce some people to confusion. That counting 
is a technical device introduced to be able to take the 
N c —> 00 limit keeping various quantities, there included 
the fermion mass, constant (unlike our result); but there 
is no reason why nature should implement this counting. 
In fact, the very concept of Grand Unification, hinted 
at by running coupling constants converging at a high 
scale, suggests that g would be the same for all groups 
(that is, independent of N c ). 

Dynamical mass generation is one of the great 
conceptual advances of the last half century, turning 
fermions that are light in the Lagrangian into heavy 
ones. The phenomenon is well known in QCD and 
we have discussed groups of larger dimension, through 
their Casimir factors Cf in the fundamental representa¬ 
tion. At the hadron scale we have employed the rainbow 
approximation of the Dyson-Schwinger equations, and 
mass generation is approximately proportional to the di¬ 
mension of the group fundamental representation, with 
a different slope for each family of classical groups (see 
figure 4). 

In the end, we have provided a plausible answer to the 
naive question Why the symmetry group of the Standard 
Model, S U(3) c X S U(2)i X t/(l)y, contains only small¬ 
dimensional subgroups? It happens that, upon equal 
conditions at a large Grand Unification scale, large- 
dimensioned groups force dynamical mass generation at 
higher scales because their coupling runs faster. Since 
the dynamically generated mass is proportional to the 
scale at which it is generated, fermions charged un¬ 
der those groups, should they exist, would appear in 
the spectrum at much higher energies than hitherto ex¬ 
plored. 


'Due, for example, to sufficiently many flavors screening the in¬ 
teraction. 
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Appendix A. Color factors 


We need two numbers from group theory, the dimen¬ 
sion of the fundamental representation of the group, N c , 
that is trivially read off, and the color factor Cf for 
fermion self-interactions, that is calculated in this ap¬ 
pendix for various Lie groups. There are two classi¬ 
cal groups for each odd N c and three classical ones for 
each even N c [39, 38]. Figure A.12 presents the result 
at a glance. It is patent that for each of the classical 


40 
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10 
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G2 
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E7 
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Nc 


50 60 70 


Figure A.12: Color factor Cf for the self-energy of a fermion in a 
gauge theory for the classical groups and the indicated exceptional 
groups. 


with C a bc the structure constants. To compute each 
group’s color factor in the fundamental representation 
Cf (its Casimir operator), we need to contract T a = 
(T a ) l j and Tb = (T b ){ from each vertex [40]. Summing 
over intermediate states (a, b, j ), 

C f 5[ - YjTa) l MT h ) ] k - J](T a yj(T a y k . (A.2) 

a,b,j a,j 

The generators are normalized by 

Tr(T a T b ) - K5 ab , (A.3) 

with k a convention-dependent constant. As we wish 
to generalize the usual SU(3) discussion to other Lie 
groups, we fix T a — with A a the Gell-Mann matrices 
and then k = 4 ■ 

The result of contracting the generators of Eq. (A.2) 
is a sum over a unique set of irreducible tensors for 
each group, either totally antisymmetric or 

totally symmetric d^ ' k , dij. k , forming a basis of the cor¬ 
responding Lie algebra 4 .We have found the following 
relations useful for the task. 


fjmf mik - aS\ , 

(A.4) 

d ijm d mjk - a8\ , 

(A.5) 


with a a normalization constant of the irreducible ten¬ 
sors, due to generalizing those to three or more indices 
(see [42]). 

Next we will study all classical groups and several ex¬ 
ceptional ones, concentrating on the very minimum and 
most important properties for this calculation and defin¬ 
ing the needed irreducible tensors. The outcome is the 
factor Cf for each group as function of the fundamental 
representation dimension N c ; the calculation is doable 
without resource to the explicit values of the generators 
and structure constants [42]. 

Appendix A.l. Classical groups 


group families, the relation between Cf and N c is linear, 
though with different slope, with the exceptional groups 
scattered and having a surprisingly small Cf for their 
large N c . 

In an iV f -colored Yang-Mills theory, the N c fermions 
q, (i = 1,..., N c ) transform in the fundamental represen¬ 
tation of an n -dimensional Lie group G, and the n gauge 
bosons A a (a = 1,..., n), in the adjoint representation. 

The T a matrices generate the associated Lie algebra 
through 

[T a ,T b ] - iC abc Tc, (A.l) 


Appendix A. 1.1. SU(N C ) 

The unitary groups are usually denoted S U ( N c ) and 
for them, 

(A.6) 


4 Note that given a tensor T e V p <8> V q , V being the vector space 
generated by a basis of p vectors, while V is its dual space generated 
by the dual basis of q forms; the tensor will have components , 
with upper indices denoting covariant, lower ones contravariant com¬ 
ponents, and both are related through complex conjugation [41]. 
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The fundamental representation of S U(N c ) is the set 
of [N c x N c ] unitary matrices with unit determinant act¬ 
ing on an /V, -dimensional complex space (N c fermions, 
for our purposes). The invariant quantities in this repre¬ 
sentation are the metric 5'j and the Levi-Civita tensor of 
N c dimensions, &i~ k . Tracing over j and / in Eq. (A.6), 



Z 


= f « 


N, 


■)»(A.7) 


and finally we reobtain the well-known result 

= (A ' 8) 

Now we repeat the calculation for other groups used less 
often in this area of particle physics. 


Appendix A. 1.2. S 0(N C ) 

For orthogonal groups SOW, ), 

- K (T a ))(T a ) l k - - 6 “6 jk ). (A.9) 

The fundamental representation of SO (N c ) is the 
set of [N c x N c ] orthogonal matrices of unit determi¬ 
nant, acting on a complex vector space which is N c - 
dimensional (for our purposes, N c fermions). 

In this representation, the invariant symmetric tensor 
is d‘i (and its inverse djj). Diagonalizing d ,j and rescal¬ 
ing the N c fermion fields q,(i - 1,..., N c ), we can always 
find a representation where d,j — 5ij. There is no dis¬ 
tinction between upper and lower indices (the fermion 
and its antiparticle), so that the representation is real. 
Tracing again over j and /, we find 

Cf = X ~(N C - 1). (A. 10) 


for N c - 2, or its multidimensional generalization. Trac¬ 
ing once again over j, l and employing the relation 

f J fjk = 4. ( A - 12 ) 

we arrive at 

C F = ^Wc+ 1). (A. 13) 

Appendix A. 2. Some exceptional groups 

Appendix A.2.1. G 2 (N c = 1) 

For the real group G 2 

- K (T a yj(T a ) l k = ] -( 6 ‘ k S l } - 6 il 6 jk ) - X -fj m r l k . (A. 14) 

The fundamental representation of G 2 (N c — 7) pre¬ 
serves the symmetric and the totally antisymmetric 
fij k tensors. It being a real group, G 2 requires no distinc¬ 
tion between covariant and contravariant indices. Trac¬ 
ing over j, 1 and applying Eq. (A.4) for the contraction 
of the fij k we obtain 

C F = i(tf c - 3) = 1 . (A. 15) 

Appendix A.2.2. E6 (N c —21) 

Next we examine the exceptional complex group E6 

+^^ 4 - l-d Um d mjk .(A.] 6 ) 

The fundamental representation of E6 (N c = 27), 
leaves invariant the totally symmetric d l ' k tensor (and 
its inverse dip). Once more, taking the trace over j, l 
and using now Eq. (A.5) to contract the d,j k tensors, we 
obtain 

1 29 13 

= -(N c -y) = -. (A. 17) 


Appendix A.1.3. S p(N c ) with even N c . 

For the symplectic groups S p(N c ) (that have the sign 
structure of Hamilton’s equations and are thus defined 
only for even N c ), 

-(TaY/Tai = hstfj - ffjk ) . (A. 11) 

K Z J 

The fundamental representation of S p(N c ) is the set 
of matrices of dimension [N c X /V„] with even N c that 
leave invariant the antisymmetric tensor f' J (and its in¬ 
verse fij), where 


Appendix A.2.3. F 4 (N c - 26) 

We now proceed to the real F 4 group, for which 

-VajjiTJt = .(A. 18) 

The fundamental representation of F 4 (with N c - 26), 
preserves the symmetric 6 ,j tensor and also the totally 
symmetric d'i k tensor. Again this is a real group, so 
covariant and contravariant indices need not be distin¬ 
guished. Taking the trace over j, l, the fundamental 
Casimir falls off in two steps. 


f U - 


0 

-1 



Cf4 = ,(A-19) 

j 


13 


c4 = ^4-4)-^. ( A - 2 °) 

C F = ^~(N C - 8 ) - 1 . (A.21) 

lo 

(This computation does require use of one explicit value 
of the totally symmetric tensor d,j k , namely that d m ^ = 
0 vanishes for a repeated index, which does not follow 
from symmetry alone). 


Appendix A.2.4. £7 (N c — 56) 
For the complex group £7, 


1 

K 


(TaYjiTai 


^(6 i / j +f il f jk --d ilmn f mj f nk UA.22) 

Z4 J Ct 


The fundamental representation of £7 (N c = 56) pre¬ 
serves the totally symmetric tensor d‘^ mn as well as the 
antisymmetric ones fy, y f IJ . Tracing the closure rela¬ 
tion over j and l. 


Cf 4 - 2 (A.23) 

j 


C ^ k =-(N c 6‘ k + 6‘ k ), 


1 57 

C F = —{N c + 1 ) = —. 
48 48 


(A.24) 
(A.25) 


(Here it has been sufficient to note that the contraction 
d l j m "f m jfnk = 0 vanishes as the tensors have opposite 
symmetry.) The color factors C F of all the groups stud¬ 
ied in this work are collected in table A.l for ease of 
reference. 


Group 

Color Factor (C F ) 

S U{N C ) 

i((V c -4-j \/jV c e N 

SO(N c ) 

|(iVc-l) VN c e N 

Sp(N c ) 

i((V r + l) N c = 2n n e N 

£6 

A (^-f) Nc = 27 

£4 

re(lV f - 8 ) N c = 26 

G2 

?(JV C - 3) N c = 7 

El 

i((V c + l N c = 56 
_ V ____ 


Table A.l: Color factors Cf for fermions in the fundamental repre¬ 
sentation needed for all the groups studied in this work. 
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